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Abstract 

We describe a scheme of constructing classical integrable models in 2 + 1-dimen- 
sional discrete space-time, based on the functional tetrahedron equation — equation 
that makes manifest the symmetries of a model in local form. We construct a 
very general "block-matrix model" together with its algebro-geometric solutions, 
study its various particular cases, and also present a remarkably simple scheme of 
quantization for one of those cases. 



Introduction 

Consider two electric devices, each made of three resistors and with three outer contacts, 
and whose diagrams are triangle and star, respectively. It is well known that they are 
equivalent iff the following relations for resistances hold: 

/ ill iii ii ?"l7"2^3 ■ 1 o o /i \ 

r,r, = r x v 2 + r 2 r 3 + r,r 1 = , j = 1,2,6, (l) 

r 1 +r 2 + r 3 w 

where resistances Tj apply to the triangle, and resistances r'- — to the star. These relations 
enable us to find r'j for any given triple of rj, and vice versa. 
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We will represent such electric devices diagrammatically in a simplified manner (in 
comparison with usual electric diagrams): resistors will be represented as edges of a 
graph, outer contacts — as blank circles, and inner contacts — as filled circles. In this way, 
the transformation ([]]) is represented in Figure 1. 




Figure 1 

Note that formulae ([3]) admit the interchange 

1 

3 

This leads to the fact that the "star-triangle" and "triangle-star" transformations, after 
a simple change of variables, can be expressed as the same transformation. The exact 



formulation of this will be given below in subsection [375 . 



Consider now the diagrams in Figure 2. We can transform diagram (a) into diagram 




Figure 2 

(b) in four steps, each being a star-triangle or triangle-star transformation for some three 
resistors. What is important, we can do it in two ways, represented in Figure 3 (where 
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Figure 3 

the obvious arrangement of blank and filled circles is not shown) . 

It can be seen from "electrical" argument that the two ways in Figure 3 must lead to 
the same result. Mathematically, this means that the transformation ([I]) is very specific: 
it satisfies (when written as in subsection |0] below) the functional tetrahedron equation 
(FTE). 

In this paper, we will try to show that FTE constitues the basis of 2 + 1-dimensional 
integrability. We will explain how to generalize the simple "electrical" construction for its 
solutions in such way that it will give not only a huge amount of new integrable classical 
systems, but even integrable quantum systems. The contents of the rest of this paper 
is as follows. In section [I] we consider the most general abstract model where the FTE 
naturally arises. In section ^] we present a concrete, but still very general incarnation 
of the abstract model, together with the algebro-geometric method of constructing its 
solutions. In section ||| we give a list of solutions representing the simplest "one-parametric 
reductions" of the general model. In section [| we show that the "electric" model described 
above is connected with the well-known bilinear Miwa equation. In section [| we present an 
amazingly direct quantization for one of our "one-parametric" models, using, essentially, 
the same FTE. Finally, we discuss our results in section |6|. 

1 Classical straight-string model and functional te- 
trahedron equation 

Let there be several straight strings, i.e. oriented straight lines (with the "positive direc- 
tion" indicated in each line), in a plane. Let there be among those strings (?) no coinciding 
ones, (ii) no three of them intersecting in one point, and (Hi) let the orientations of strings 
be consistent in the following sense: moving from any point along the strings in positive 
directions, one cannot return in the same point. Let us attribute to each point of intersec- 
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tion of two strings, e.g. strings a and b, some object X ab . Now let us allow each string to 
move, remaining parallel to its initial position, and pass through points of intersection of 
other strings, but so that (i) two parallel strings never go through each other and (ii) at 
no moment there should be four strings intersecting in one point. 

Let us assume that an object X ab changes only if some string c goes through the 
intersection point of a and b. To be exact, in such case the triple (X ab , X ac , X bc ) is 
transformed into a new triple (X' ab , X' ac , X' bc ) according to some fixed rule. When we 
consider such a transformation, we think of both those triples as ordered in the following 
way: one can get from X ab to X ac and then to X bc moving along the strings in positive 
directions, while from X' ab to X' ac and then to X' bc — in negative directions. This can be 
represented graphically as turning the triangle inside out, see Figure 4. 




Figure 4 

We will say that the operator R acts on the triples of objects: 

R : (X ab , X ac , X bc ) — >• (X' ab ,X' ac ,X bc ). (2) 

Note that in Figure 4 the point of intersection of strings a and b goes from the left part 
of the plane, with respect to orientation of c, to the right part. If we make in Figure 4 
the changes c <-> a, X <-> X' , LHS <-> RHS, we will see that for the case when the point 
of intersection of strings a and b goes from the right part of the plane, with respect to 
orientation of c, to the left part, R should be replaced by P^R^P^, where P± 3 stands for 
the interchange of the first and the third objects in a triple. In other respects, we assume 
R be the same for all triples of strings, i.e. not depending, say, on the angles between 
them or whatever. At the same time, we will not exclude the situation where R is not 
single-valued — in such case we will say that the "primed triple" is determined up to a 
"gauge equivalence" . 
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If there are many enough non-parallel strings in our model, there exist many possibil- 
ities for their movement. It is natural to regard as "integrable" such a model where to a 
passage from one string configuration to another there will correspond a transformation of 
objects X a b not depending of the details of that passage, e.g. which string was the first to 
pass through the intersection point of two others and so on. For this to hold, it is enough 
to require that the operation of turning a triangle inside out, as in Figure 4, commute 
with passing a fourth string across the whole triangle. In other words, the two ways of 
transforming the left side of Figure 5 into its right side, one of which starts with turning 




Figure 5 

inside out triangle 356 and the other one — triangle 123, must lead to equal results. This 
condition is the functional tetrahedron equation (FTE) to which this paper is devoted: 

#123 ° #145 ° #246 ° #356 = #356 ° #246 ° #145 ° #123- (3) 

Here, of course, #123 acts only on the triple 

(Xi,X 2 ,X 3 ) = (X ab , X aC} X bc ) 

and so on. 

The diagrammatic representation of equation (||) itself is given in Figure 6. 

To conclude this section, note that we could assosiate X-like objects not to the inter- 
section points of strings, but to the segments in which the strings divide each other, or to 
the domains of the plane in which it is divided by the strings. It would lead us to other 
versions of the functional tetrahedron equation, in complete analogy with the quantum 
case [pi 0]. It is natural to call equation (^) the vertex type FTE. In this paper we will be 
dealing only with this type of equations. 
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Figure 6 

The block-matrix model and its algebro-geometric 
solutions 



2.1 Formulation of the model 

Let now to each of the strings correspond a finite-dimensional complex linear space, e.g. 
to the string a — the space V a . We will also associate with a set of strings the direct sum 
of their linear spaces (and not the tensor product, in contrast with the usual theory of 
Yang-Baxter equation). The objects X ao will be linear operators in V a © V&. Fixing once 
and for ever the bases in all V a , we will identify those operators with block matrices of 
the type 

A B 
C D 



X 



where A acts within V a , B- 
the action of operator R (||) 



-from Vb to V a , etc. The "primed" operators- 
— will be determined from the equation 



X„hX nr X\ 



ac-^-bc 



Xb c X ac X ab , 



-the result of 



(4) 



which in expanded form is written as 




/l 

A' 3 B' 3 
VO C> 3 D' 3 



(A> 2 B> 2 
10 

\C 2 D' 2 



10 

A 3 B 3 
C 3 D 3 

A' x B[ Q 
C[ D[ 
1 



(5) 
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i.e. we identify each of the operators 



Y - Y - ( Al Bl \ Y< - Y' - f A 3 B 3 

^"^"id Dj ' ' ' ' ' Afoc ~ A3 ~ \C 3 D' 3 

with the direct sum of itself and the unity operator in the lacking space. 

Remark. If we now want to consider Figure 4 as the illustration to formula @j, we 
have to take into account that the operators in (H) act against the arrows in Figure 4. On 
the other hand, in Figure 6 the arrows point in the same direction where the (nonlinear) 
operators R act. 

In case if the operators in the LHS of equation (^) are generic enough, that equation is 
solvable with respect to the "primed" operators, and those latter are determined uniquely 
up to the following obvious block-diagonal "gauge transformations" : 

A> 3 B' 3 \ (A> 3 B' 3 \fL 



c 3 d' 3 j y&t d' 3 j \ o m 

A' 2 B' 2 \ _ f 1 \(A' 2 B' 2 \(K 

C 2 D' 2 ) ~* \0 AT 1 AC£ D' 2 )\0 1 

A[ B[\ (K- 1 \ (A[ B[ 

C[ D'J { L"V [C[ D[ 

The details can be found in 0], chapter 2. 

As for the functional tetrahedron equation, it follows in our block-matrix model from 
the fact that the product of six X-matrices corresponding to the RHS of Figure 5 deter- 
mines those matrices as well (again up to obvious block-diagonal gauge transformations). 
The proof of this statement is a routine exercise in matrix algebra (to begin, one can 
unite, for a while, the strings c and d in the RHS of Figure 5 in one "thick" string, and 
then apply several times the uniqueness, up to gauge transformations, of factorizing the 
product of three matrices). 



2.2 Algebro-geometric solutions 

Let us take some (smooth irreducible) algebraic curve V of genus g. Let us attach to each 
string a two effective divisors (i.e. finite sets of points of T with positive multiplicities) 
T>' a and D" a of the same degree, equal to dim V a (the dimension of linear space associated 
with a). For two non-parallel strings a and b, we will require that no point enter at the 
same time in both T>' a and T>'l (as well as in both T>" a and T>' b ). The meaning of this will 
be clear later. 

Next, let us put in correspondence to each point A of the plane where the strings live a 
divisor T> A according to the following rules. First, for each string a, if A lies to the left of 
a (with respect to its orientation), take divisor V' a , if A lies to the right of a, take divisor 
V", if A lies in a, take the zero divisor, and denote that as Then set V A = J^aii a^a- 

We will only be interested in divisors corresponding this way to (inner points of) string 
segments limited by string intersection points, and to those intersection points themselves. 
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For example, if there are only two strings, to their halves and to the intersection point 
will correspond divisors as in Figure 7. 




Figure 7 

Take now an arbitrary (generic) divisor T> of degree 

deg V — ^2 dim V a + g — 1 . 

over all a 

If to a given segment / of a string a corresponds a divisor T>i, let us also put in correspon- 
dence to that segment the linear space of meromorphic functions / on F whose divisor 
(/) (the difference between the zero divisor and the pole divisor) obeys the condition 

(f)+V-V l >0 

(this means that the singularities of / are controlled by divisor T> and, moreover, / must 
have zeroes in the points of divisor Dj). It is easy to see that the space of such functions 
is isomorphic to (i.e. has the same dimension as) V a , so we will identify it with V a , using 
for that an arbitrary isomorphism. 

We will need the following explicit description of isomorphisms of this kind. The 
space V a itself, whose basis, as we have agreed in the beginning of subsection [2.1| , is 
fixed, can be imagined as comprised of column vectors (that we will write sometimes as 
transposed row vector, using the symbol T of matrix transposing) of height n a = dim V a . 
Whenever V a should be identified with some space of meromorphic functions, we should 
choose in the latter a basis (fx, . . . , f na ), and after that put in correspondence to a column 
vector (ai,..., a na ) T <E V a the function a 1 f 1 + ... + a n J na . 

Thus, in the different segments of string a, the space V a is identified with different 
spaces of meromorphic functions by arbitrary but fixed isomorphisms. Every time as, 
following the movement of the strings, the limiting points of a segment pass through each 
other, the isomorphism for the "turned inside out" segment is chosen, of course, anew. 
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Let strings a and b intersect in a vertex A. Let us denote Da the divisor corresponding 
to A according to the rules described above. The vertex A has two incoming edges (string 
segments preceding A, according to the given positive directions in the strings) and two 
outgoing ones. It is easily seen that the sum of linear spaces of meromorphic functions 
corresponding to the incoming edges consists of such functions / that 

(f)+V-V A >0 

and is direct (the particular case shown in Figure 7 can be useful to quickly understand 
this. Recall the requirement that no point enter in both T>" a and T>' b ). Exactly the same 
applies to the outgoing edges, hence the sums of "incoming" and "outgoing" spaces coin- 
cide. 

The operator X a b is now constructed as follows. Take some vector $ G V a © VJ, and 
represent it as a sum 

$ = $ a + $ 6 , $ a G V a , $6 e V b . 

Identify vectors $ a and with meromorphic functions (p a and (pb using the discussed 
above isomorphisms for the outgoing edges. Thus, vector $ is identified with function 

<f = + <fb- 

Now decompose this function as 

where ip a and ipf, belong to the spaces corresponding to incoming edges, identify ip a and 
■06 with vectors ty a G V a and $(, e 14 respectively using the isomorphisms for incoming 
edges, and set 

= ^ a + 

The operator X a b is defined by the requirement that 

X ab § = * 

for every vector $ (recall that X ab acts "against the arrows"). 

The equation ([|) for X-operators constructed this way follows from the fact that both 
LHS and RHS of (§J) act on each "outgoing" vector $, where now $ G V a @V h ® V c , 
as follows. First, $ is identified with a meromorphic function ip using isomorphisms 
corresponding to three outgoing edges. Then ip is identified with some vector $6^9) 
Vb ffi V c using three "incoming" isomorphisms. So, both LHS and RHS transform $ into 

although using different ways of doing this. 

The explicit formulae determining X ab are as follows. Recall that the isomorphisms 
between the spaces V a of column vectors and the spaces of meromorphic functions are 
provided by fixing bases in the latter spaces. Let (fx, ... , f na ) and (#!,..., g n J be bases 
corresponding to incoming edges for some X a b, while (fi, ■ ■ ■ , fn a ) an d Q?i, ■ ■ ■ , <7nJ be 
bases corresponding to the outgoing edges. Then it is not hard to verify that X ab is 
detrmined by the formula 

(h(z), • • • , fn a (z),gi(z), . . . , g nb (z)) X ab = (fi(z), fn a (z),gi(z), . . • , g nb (z)) 
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that must hold for all z G Y. 

Remark. Our model, being composed of a finite number of strings moving in the 
plane, has the obvious "global" integral of motion (or rather all admissible motions) — the 
product of all X-operators (complemented by unities in the lacking spaces) taken in the 
order corresponding to moving along the strings in positive directions. If, however, we 
want to ensure the possibility of passing to the infinite string number limit, we should 
prefer more "local" constructions. For example, using the presented method of construct- 
ing algebro-geometric solutions, one can add a new string — let us denote it h — so that it 
don't intersect some chosen domain Q of the plane, adding at the same time to divisor 
T> one of the divisors T>' h or T>'^, depending on whether Q lies to the left or to the right 
of h. The result will be no change to the X-matrices (as well as divisors corresponding to 
string segments) in Q. The reader can find some formulae for the infinite kagome lattice 
in [|J, chapter 3. 



3 One- parametric solutions to the FTE 

The simplest case of the "re-factorizing equation" (|5|) is when the blocks in matrices are 
just numbers. For numbers, we will use small letters, i.e. let 
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b' 3 
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and let us search for solutions of equation 

^12^13X23 = -^23^13^12- (6) 

The (nontrivial) one-parametric solutions of (||) were classified in 0. The word 'one- 
parametric' means here that a^, bk, Ck and dk, as well as a' k , b' k , c' k and d' k , are supposed 
to be algebraic functions of one complex parameter x in such way that 

a k = a(x k ), b k = b(x k ), c k = c(x k ), d k = d(x k ), 

a k = a(x k ), b k = b(x k ), c k = c(x k ), d k = d(x k ) } k = 1,2,3, 

and the operator R can be regarded as a transformation 

R : (x 1 ,x 2 , x 3 ) — > ( 
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It turns out that there are, essentially, six different cases for the matrix 

X(x)- 

3.1 Case (a) 



a(x) b(x) 
c(x) d(x) 



k being a constant; this gives 

R : xi, x 2 , x 3 — > Xi, kx 2 + x x x 3 , x 3 . 

Inverse map: 

D _l X 2 - X!X 3 

3.2 Case (/?) 

X(x) 

this gives 



1 a; ' 

fc/x o ; ' 



fcc 2 + Xix 3 kx 2 x 3 
R : Xi,x 2 ,x 3 ^> , x ± x 3 , 



x 3 ftx 2 + £i£ 3 

Inverse map: 

R' 1 : xi,X2,x 3 



x 1 x 2 X1X3 x 2 + Xia; 3 



x 2 + X1X3 A; xi 
3.3 Case (7) 

*<*> = (if, ?)• 

then 

x 3 - x 2 + xia; 2 xix 2 a; 3 

i? : xi,x 2 ,x 3 ^ , ■ , x 3 . 

x 3 x 3 - X 2 + XiX 2 

Inverse map: 

D-l X 1 X 2 

R : X 1 ,X 2 ,X 3 ^ ; , X 3 - XiX 3 + XiX 2 , x 3 . 

X 3 - XiX 3 + XiX 2 
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3.4 Case (5) 

X(x) 

Then 



x 1 
1 -x 



XiX 2 (1 — 2])2 2 2 3 

i? : 21,22,23 -> ; , 2i + 2 3 -2i2 3 , 



2i + 2 3 - 2i2 3 2i + 2 3 - 2i2 2 - 2i2 3 

Here 

R 2 = 1. 

This transformation is connected with the pentagon equation and was described in 



3.5 Case (e) 

X{x) 

then 



2 1 + ix 
1 — ix x 



R: X 1 ,X 2 ,X 3 —> ; ; , 2i + 2 3 + 2i2 2 2 3 , ; ; , (7) 

2i + 2 3 + 2i2 2 2 3 2i + 2 3 + 2i2 2 2 3 



again with 



R 2 = 1. 



This is the electric network transformation, considered in |]J . The variables Xj and x'j are 
connected with the resistances r 3 - and r'- in the Introduction as follows: 

1 1 , , , 1 , , 

X\ — — , 2 2 — r 2 , 2 3 — — , x 1 — r 1 , 2 2 — — , 2 3 — r 3 . 

n ^3 ^2 

The reader can verify that the equalness of results of two transformation sequences in 
Figure 3 means exactly the same that the FTE for x's, if all the resistances attached to 
the edges of graphs in Figure 3 are replaced by xf l , where the proper choice of signs is 
an easy exercise. 

Note also that in |L| this transformation was realized in terms of a "free-bosonic local 
Yang-Baxter equation" , while our case here is a fermionic one, as explained in 0. Q . 



3.6 Case (C) 

where s(2) 2 = 1 — x 2 . This case is equivalent to Onsager star-triangle transform, and can 
also be interpreted as the Euler decomposition of an element of the group £0(3) into a 
product of two-dimensional rotations, if we put x = cos0, s(x) = sin (p. This case was 
also considered in [jlj. 
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For R, this gives 



i? : £1,22,23 7,^(^1,^2,^3; 



F(2i,2 2 ,2 3 ) F(Xi,X 2 ,X 3 ) 

where F can be found from 

s(F) = s(x 2 )x 1 x 3 - s(xi)s(x 3 ), 

and so this is a two-foiled transformation. 
Inverse map: 

R- 1 = IRI, 

where I transforms a pair (x,s(x)) into (x, —s(x)). 

4 Electric network transformation and Miwa model 

It has been shown in paper Q that the "electric" model described in the Introduction and 
in subsection |375| is connected with the well known integrable Miwa model H]. In terms 



of the present paper, this connection reads as follows. 

Let us consider, in the three-dimensional space R 3 3 (x,y,z), the cubic lattice planes 
whose equations are x = I, y = m, and z — n, where l,m,n take all integer values. 
Consider also a "moving" plane with the equation 

x + y + z = t, (8) 

where t can be thought of as "time" . With a generic t, the intersection of the cubic lattice 
planes with the plane (|^) yields in this latter the regular (infinite) kagome lattice. Let us 
choose some positive direction for every straight line forming that lattice in such manner 
that all parallel lines be directed the same way and the orientation of all the lines be 
consistent in the sense of section [l| 

Let us now attach to each kagome lattice vertex a variable x/~, and require that the 
strings move with time t according to equation @ and the variables Xk change according 
to (0) whenever a vertex passes through a line. In such way, we will obtain a straight-string 
model of the kind described in section [I]. 

Let us introduce unit vectors fx, f 2 , f 3 pointing in the directions of axes x,y,z of 
the space R 3 respectively. Note that each kagome lattice vertex sweeps between its two 
collisions with strings exactly an edge of the integer cubic lattice in R 3 . Thus, we can 
consider values x^ as corresponding to edges of the cubic lattice. For convenience, and 
also in order to link our presentation to the paper [fij, we will think of x^s also as attached 
to vertices of the cubic lattice as follows: let the radius vector of some vertex be n, then 
we will denote as Xi(n) the value x\ corresponding to the edge linking the vertices with 
radius vectors n and n — fx, as x 2 (n) — the value x 2 corresponding to the edge linking the 
vertices n + f 2 and n, and as x 3 (n) — the value x 3 corresponding to the edge linking the 
vertices n and n — f 3 . 
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Let us attach signs e(n) = ±1 to the cubic lattice vertices in a checkerboard order. 
Then let us introduce (again in order to link our paper to the work |]) the vectors e\, e 2 , 
according to 

fz — f'2 — fl fl + f2 + fz fl — J2 — f?, / n s 

ei = 2 ' 62 = 2 ' 63 = 2 • (9) 

Finally, fix some non-zero constants ai,a 2 and and make a substitution 

*,(») = «(») a, r(« + e .+ e , + «.)r(n - ej ) 
r(n + e fe )r(n + e;) 

where r is the new unknown function, defined, as is seen from in the centers of cubes 
of the lattice. 

As is shown in U, function r satisfies the Miwa equation 

4 

ctjr(n + ej)r(n — ej) = 0, e4 = — ei — e2 — e3, «4 = ai«2«3- 

Note that the variables Oj(n) in |j] are connected with our x 3 -(n) by 

% (n) = e(n)a; i (n)~ 1 . 



5 Quantization of the (/?) model 



The model corresponding to the case (J3) (subsection |3.2|) admits an amazingly straight- 
forward quantization. According to subsection |3.2| , the corresponding transformation 
(xi,x 2 ,x 3 ) — > (x[, x' 2 , x' 3 ) arises from the relation 

1 x 2 
1 

k/x 2 



/ 1 x l 


0> 


k/xi 





V o o 




/l 




, 1 


i) 


\0 k/x' 3 





1 x' 2 
1 
k/x' 2 




(10) 



To quantize this model, let us take the same relation (|T0|), but with Xi,x 2 and x 3 
belonging to some associative algebra and satisfying the commutation relations 



XiX 2 = lo l x 2 x\. 



X\Xz = LO X3X1, 



x 2 x 3 = u) l x 3 x 2 , 



111 



where u is a scalar (and k remains a scalar, too). It can be verified that in this case (JX0| 
still determines the "primed" variables without contradiction, namely 



x\ = x\ + kx 2 x 3 1 , 



.X n 



+ kx 2 x 3 ) 



x 2 , 



(12) 
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and the following relations hold: 

*^1*^2 — ^ ^2*^1 5 X-^X^ — UJ X^X-^^ X^X^ — UJ X^X^. ( ^~^) 

Let us consider now a sraight-string model, as in section [I], and define the commutation 
relations for Xk in all vertices. We will demand that 

• if two Xk don't lie in the same string, they commute; 

• if the vertices where two Xfc £1X6 situated belong to a triangle of the form as in 
Figure 8 (a) or (b), the commutation relations (|ll]) hold (where we temporarily 




Figure 8 

change the actual /c's to 1, 2 and 3). Note that there may exist other strings, not 
depicted in Figure 8, between the vertices 1, 2 and 3; 

• if the vertices where two x^ are situated arise as intersection points of a string with 
two parallel strings, the commutation relation between them can be obtained if we 
let one of the vertices in Figure 8 tend to infinity. Here the figures (a) and (b) lead 
to the same result, however the orientation of strings does play the role. 

It can be verified that the stated commutation relations are conserved under the 
transformations as in Figure 4. 

6 Discussion 

In this paper, we have considered, first of all, classical integrable dynamical systems in 
2 + 1-dimensional discrete space-time. We proposed a very general scheme for constructing 
such systems together with their solutions. We showed that our scheme includes Miwa 
bilinear equation as just one of particular cases. 
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The integrability in our general scheme is ensured by a rich set of symmetries that 
are expressed in the local form as functional tetrahedron equation (FTE). The form of 
FTE suggests at once that the corresponding model need not be associated with a regular 
(e.g. cubic) lattice, but can be naturally defined e.g. for any configuration of intersecting 
planes ("world sheets" of strings in section |I|) in a three-dimensional space. The algebro- 
geometric solution for such system is given in terms of an algebraic curve and some 
constant divisors assigned to the planes. In a sense, the equations of motion, when 
expressed in terms of those divisors, "dissappear" , and this calls to mind the Penrose 
twistor theory, where field equations "dissappear" as well. 

We have also shown that, at least in one particular case, FTE can successfully replace 
the usual (quantum) tetrahedron equation in constructing a quantum integrable model. 

Finally, we would like to note that a generalization onto 3 + 1 dimensions of our 
scheme for constructing classical models and their algebro-geometric solutions has been 
found recently in work [^]. 
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